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SUMMARY 

The Batdorf model is modified to include the reduction in shear due to the 
effect of compressive stresses on the crack face. This new formulation was 
used to obtain the probability of failure of ceramic components under contact 
stress conditions. The combined effect of the surface and volume flaws are 
included in the analysis. Due to the nature of the fracture of brittle materi- 
als under compressive loading, the component is modeled as a series system in 
order to establish bounds on the probability of failure. 


INTRODUCTION 

The problem of failure of two bodies in contact goes back to the period 
when railroad wheels, gears, and bearings were first used. The classical 
deterministic approach to design against failure is well accepted and widely 
used for metal components. However, because of the attractive physical and 
mechanical properties of modern ceramics (lightweight, high-temperature 
strength, and wear resistance), ceramic components have in recent years been 
considered for structural applications. As in the case of other brittle mate- 
rials, ceramics exhibit a large variation in fracture stress that must be taken 
into account in design. This variation in strength results from the presence 
of microscopic random imperfections or flaws. Ceramic components contain two 
types of flaws: volume flaws and surface flaws. Volume flaws arise from mate- 

rial processing while surface flaws arise from grinding and other surface fin- 
ishing operations. 

Most of the probabilistic approaches to brittle fracture are formulated 
for tensile failure. However, ceramic components are being used under compres- 
sive loading, as in bearing applications. In this paper, the probability of 
failure under compressive loading is formulated based on the Batdorf tensile 
model . 


REVIEW OF THE BATDORF MODEL 

The failure of brittle materials is attributed to the presence of micro- 
scopic flaws. The material fails when the strength of the weakest flaw is 
exceeded. The weakest link model was proposed as a statistical theory of frac- 
ture (ref. 10). These flaws were assumed to be cracks whose strength was 
dependent on their size and orientation (ref. 3). 

Batdorf assumed that these cracks were uniformly distributed and randomly 

oriented as shown in figure 1. The materials will fail when the effective 

stress on a crack reaches some critical value cr characteristic of that 
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particular crack. To determine the effective stress, assumptions are made con- 
cerning the shape and the fracture criterion (ref. 4). 

Two of the mixed-mode fracture criteria are the maximum tensile strength 
and the strain-energy release rate. For this analysis, the strain-energy 
release rate criterion was selected because of its greater degree of shear- 
sensitivity. The effective stress applied on a crack is given as 


for a Griffith crack, where is the macroscopic stress normal to the plane 

of the crack and t is the shear parallel to it. A crack will fracture when 
the effective stress is greater than or equal to the critical stress. 


COMPRESSIVE MODEL 

When the stress normal to the crack plane is compressive, the shear and 
the friction due to the normal stresses act against one another (ref. 1) 

0 = T = T + ua (2) 

e e n 

where x g is the effective shear on the crack and y is the internal friction 
coefficient of the material. 

To determine the probability of failure, the material is divided into ele- 
ments. Within each element the stress state is assumed constant. Each element 
contains a number of cracks. The crack density N is material dependent and 
is defined as the number of cracks per unit volume having a critical stress 
greater than or equal to <7 cr . This function is determined by experiment using 

uniaxial data (ref. 6). The probability of failure for volume flaws is 

P fV - ’ - P SV ■ ' - **<> (-!„ dv f" m " h d °cr) <3 > 

where n is a solid angle containing all of the orientations for which 
a > a cr . Then the quantity Q/4ir is the probability that the crack will lie 

within the angle since all orientations are considered equally likely. An 
example of this angle is shown in figure 2. A similar approach is taken to 
find the probability of failure for surface cracks (ref. 5) 



where N s is the surface crack density function, A is the surface area, and 
a is the radian measure of the angular range in the positive 0 ^ half-plane 

within which o fi > o cr . The total probability of survival, P$, is equal to 
the probability of survival for the volume times the surface survival, 
(P SV P SA)- 


2 



The upper limit for the integration of o cr is the maximum effective 
stress CT emax - When CT cr > CT emax’ ^ or u equal to zero. For tensile 
stress states the maximum effective stress is equal to the maximum principal 
stress. For compressive stress states, the maximum effective stress is depend- 
ent on the stress state and the internal friction coefficient. If o < 0, the 

e 

frictional force along the crack is greater than the shear and the crack is 
locked. This phenomenon is taking place near the orientation of minimum prin- 
cipal stress as illustrated in figure 2. 


SYSTEM RELIABILITY UNDER COMPRESSIVE LOADING 

Brittle materials with pre-existing cracks may fracture when loaded in 
triaxial compression (ref. 2). The unequal principal compressive stresses gen- 
erate shear stresses which act against frictional forces, initiating local 
crack growth. Since the existing cracks are microscopic, a single crack does 
not produce total failure as it almost always does in tension. Total fracture 
occurs in compression when several of these cracks extend and join together 
(single crack arrest is a possibility under high compressive loading). To 
account for the multicracking phenomenon, the material is modeled as a series 
system rather than independent elements as in the weakest link theory. Limits 
are established on the probability of failure of the series system using the 
Ditlevsen bounds (ref. 7) 
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(5) 


where P^ is the probability of failure of an individual element, assembled in 
decreasing order, k is the number of elements and Pjj is the joint probabil- 
ity of failure of the elements i and j, which is formulated for simplicity 
in function of the safety indices |3. and |3j , given by 


p ij- p iVJ'o' p< -V-V z> 


dz 


( 6 ) 


where 


<p(x,y:p> = ; exp 

( 2 „V7T7) 

The safety index (3. = -4>~ ^ ( P ^ ) is a normally distributed function and p is 

the correlation coefficient (p = 0 no correlation and p = 1 fully correlated). 
When p is equal to one, a k series system is modeled as one single element 
whose probability of failure is the average of the k elements. 


1 (x 2 ± y 2 - 2xyp) 


(1 


p 2 ) 


(7) 
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APPLICATION 


The volume crack density function, N in equation (3), is assumed to be 
(ref. 3) 


where kp and m are material constants. A similar expression, N„ = k D a m , 
w s ds cr 

was used as the surface crack density function in equation (4). These crack 

density functions were calculated based on compressive experimental data of 

alumina (ref. 9). To determine the crack density function the probability of 

failure versus the absolute value of the compressive strength was plotted. The 

resulting curve may be expressed (ref. 10 ) 


P 


f 


= 1 



( 8 ) 


where o is the compressive strength and C and m are constants determined 
by linear regression. The constants kg or k Bs are determined by equating 
the exponent in equation ( 8 ) to the one in equation (3) or (4), respectively, 
for example 
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where m = 43.41, k B = 1.482xl0 5 GPa (m)0-069 and kgs = 3.177 GPa (m) 0.046. 


The preceding theory was used to evaluate the probability of failure for 
an alumina ceramic under contact stress loading as shown in figure 3. A closed 
form solution to determine the stress distribution beneath the contact region 
was used (ref. 8 ). The probability of failure as a function of the maximum 
contact pressure is shown in figure 4 with p = 0.35. 


When the elements are not correlated the weakest link probability lies 
within the Ditlevsen bounds. However, when the system is fully correlated, 
the probability of failure is substantially lower than the weakest link model, 
for a given load. This indicates that for the whole system to fail under com- 
pressive loading, more than one element needs to fail. An intermediate corre- 
lation function was assumed p.^ = exp -(r^j/d - 1 ) where r^j is the distance 

between the centroids of two elements i and 3 and d is the mesh size. 

This correlation is more realistic because p is equal to one for two adjacent 
elements and decreases as the distance between elements increases. 
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FIGURE 2. - SOLID ANGLE WITHIN WHICH THE NORMAL TO 
THE CRACK MUST BE ORIENTED FOR FRACTURE TO OCCUR. 
( 0 1 = 0 2 = 0 , 0 3 < 0 .) 



FIGURE 3. - SCHEMATIC VIEW OF A CONTACT STRESS DISTRIBUTION ON A 
SEMI-INFINITE REGION. 



FIGURE - PROBABILITY BOUNDS OF FAILURE AS A FUNCTION 
OF THE NORMALIZED MAXIMUM PRESSURE FOR DIFFERENT 
CORRELATION COEFFICIENTS, (p). 
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